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SYMBOLIC IMTEGRATION 
by Joel Woses 

t * '* ^f^^^^lT^ ^®^ ^^^ written which Is capable of 
Intjuratlrif all hut tm of the problems solved by Slaele's 
symbolic Intmtratlon proirann SAIHT. In contrast to SAINT, It 
Is a purely alprtthiilc program,, and It has achlaved rmnlnl 
t mes two to thr«€ ordart of macultwda faster than SA?Nt! 
Tl I s program and spm« of tha basic routines which It uses 
are da.scrlbad. A haiirlitic for Intntratlon, called the Fd« 
*'r[^!J'''^i!* presented. It Is claimed that this heurlstr 
Mith the aid of a few alforithmf Is capable of s^lvlZ aH 
otL^ai^^Zo '^* ■tlforlthmlc program 1,7 rr^.; 



Sla«le*s program SAINT (Symbolic Automatic iffTegrator) 
which was flnishecl In 1961 Is cttrrently the onlf piiblisfied 
general Integration prograra, !t !s unquestlofiabiy a tour de 
force In recursive progranmil ng^ In the use of heuristics , 
and 111 the simulation of hitman behavior by a compiiter. 
However^ many people wlio have read the description of S/UfIT 

have obtained the Irapresslon that one raijst use hem istic ii Ii j 

order to solve Integration problems on the level of SAIM r 

C !*€:■, fresiifiian calcultis eitam pro'blems) „ Me shall corr'-tlder 

tiere two approisclies to Integration, Flr$t,# we find that h ^ 

order to^ perforrn Inttgration problems on the level of SAINT 

i^ni even better)j, one can ' do very nicely with a i:)i li 'ely 

algorfthnilc pr-agraai. Second,, we consider haorlstlcs 'which 
are necessary in order to be able tO' haindle those cp^<^f-% 
which fall outside thcs range of the algorf tliiiis. 

Me haye exper I merited with a proigrarn which has ninr^ 

algorithms and have^ found It to be able to solve all but t-w i> 

of the problems that SAIMT solved. It was able to solve" nnp 

of the ti«_ problems atteropted but not solved by SAIMT. ! i: It , 

also capafcile of solving many problerns not attempted" hy 
SAIMT. We do not believe that the converse Is true to nnv 

appreciable extent. The program runs at speeds m hich -ii p 

freciuently two. to threes orders of magnitude faster than 

SAIWT even thougfii nwch of the pi ogram Is still yrisomprii?. 1 

and chaining takes half the running tlrw*, lie havft ' also 

done some hand sliiBjliatlon on a lieurtstlc for iritepratlo !■ 

which appears remarliiabl; ?■ powerful* The heuristic Is « «- ijo<v=; 

at thcij t'orro of thi:! solution^ differentiate the foniL, ' aad 

deternriirii;! ths v. lines of the variables In the form (ms, !;!!■■ ly 

^¥ m '::ssl! ig sgalnl. The hetirlstlc (called the Edge hpurl v^ir 

- for iiOticated QiiEss) Is an eKtenslor* of the" rrifithfv:: f )f 

I ntegrat Ion by parts. it Is recursl ve>. and ten* Is i o 

gciiicrate si ii3i:]Too'ieiii:s which an ! progress I Visly simpliprr, l-Ji'i h 

this I'tetiristlc^ f<i/ irnifliici, we can solve the problams ti sat- 

our algorithfiilc -ir ^j iMi.t uas Uiiiabla to solve and which f ii; ili\| r 

i=ifasc:, it Is oyr n^jpo iliai^-iftfi heitrl sties such as *thi'i''"vf« 

^^'^^^ ''' • 3i^>'ie *;«' 1/'^' tin solution to problems wl ilrh" ari« 

tr«ily m ff Iciiif i\ 1, u fioriJtlcSor Intograls In sts^ idard 

!njo«t'ai;lori tablifM' u Up f-ilrly cT'mple and obtali labla u f r h 

'f*:'»?P'»''2 ' '-ase v?'tli i auuill number of algorf thras which ^ri» 

C£'^;ilv codedc Tir; j t I'orl clujis for the lritei?rratIoni pmrra-^j 
>'C)7ft fjiciecii m4 d<5tNp; '(•, i In dboiit two wseksc 



"1 " h a 3 1 go r I Z. h m s w h ! c h h a v e l;> e g ri 
bc'iOi*!,, Here R(z,) arid SCz) stand fo 
(ratio of tvio pa I y nom I a Is) ! n z ,. 

I) a) x''^""' R(;c" } syi;jsl:!ttite !/- k '^^ 

J K ''"/'Ck''^ -MDcf] [ tjeccjmes/y/Cy-^' «■!) !y 



i!::?,) 



b) IUk) ciialns to a proip'aiTi wri t:!:©!! by M'^}^. ms e'" v:"r: ■ ':■ 

Integrates rational functions by cj meiir't :)d die:, ■;: rlbed i^ 

Hardy. 

I i ) R C eiifii X ) ) ;■:: y!>s: 1 1 1 ti te y«exp ( ;{ ) 

j e K p C X ) / C 2 « 3 e k d ( 2 ;i ) ) <l ;; t:i ecom e s / 1 / ( 2 •■:- 3 y '^ !• d ': ;" 

lit) RCx ^ ya-^bx) !! ;5jbst:it:tJt€? y «/a -Msx 

/ K /x*> 1 clx becciTias yS: ( / • I } /d y 

T h i s p rc) i> I « »m con 1 cl n o £ be s o I ve g b y 3 A I lU ,. 

IV) RCk ,iB +• t«t') 

I ) a /"O .p b :> (I s i i b s I: i t u t e t a i i C y ) = x 

/ 1. / C h * >i: ■■■'' ) d K t } ecmm s /l / 2 s e s { y ) d y 

1 1 } a >CI ,, b-< ;= , I I b s t f t u t c s I r-^ C :* ' ) « ; .; 

/ K / C I "■■ ic "' ) ^ cl ; ■ b e ci mte s fs I f; C y ) d ! ?■ 

1 1 ! I a < ,, b 7 3 Ui b s 1: 1 1 u t & :: , e c C y ) =' y 

I ) c '"" s II ,b s 1 1 1 1.S t e y » jc k -r- lA:: ?: ' i* b .?.: *■ 3 

1 1 i' a ,? s ( I b s t ! 1: « I: e ( / c »'■ • "'■■■ b k < a • ■ 7a ) / n 



V i ) i-i C s I F C *'■: C H ) ) ,, « ih: r s '■* e 1 1"; g F ;;:'■; i: c ,,;•: r 1 1; o -f !' ■;: C x : I 

also* 



^at'b ;:•!;■':! 



1 1 F=- log 

S o "I u t: i ■ ■) ! ■; I s IF ( S C y ) ) / H i x ) cf ^: • j k C y ) cl k !> t ^ } / S ( ^k : :,:: ;k, 

yx i og ( y ) d 2C i:5ec:on]s s 1 / 2 x;" ' 1 c?g C x } « / ;t / 2 d s 

I ! ) ¥■•- Brctfsn 
So 1 tit 1 1:5! r I 3 F C 3 C y } } j R C y ) dj<: -;/iR C k ) diyS { y ) / C I ■■?■■ S'''' C u ) ) ..:: ::c 



'i '"■: a i^ <■■• ''3if*CS j Ti 



So 1 tit I '!:,;» 51 I :!:; f { $. f v; "!:; ) I ! I ( Tf '^ ,: ? J ,.. 

■ "" 1 • • ' • '■■■- ■' ■■ 

I 7 , ..J ,, 

X a rca I r ( .k I'd ;-i becci-!Tiie s I / 3 .k " ' e r i;:;, 

Vi I ) , "l!:CF; (a ■<'•!??!:,)) v^:,,?! ;, hi-: ; i :h 

P r €i p a 6" .3 1 o r y s t e p : y » c -< • l:i i; y ! .:3 1 



J i i! >;; I' !'• 



i ) T f s n £5 f X3i t"! Ti a 1 1 t: ir I go 5 i Ofnif 5 1 ? '* ■;":: f ' i * i ;; 

cos In©::'; and Zest to see If sIk? r^ ^■■; m r:": ' | 

■| ^,„,. , ,,;«■' -i,, __^^, , ,^^ __3 _, 



s 1 11 C f ) CI y b e c: 'srssiiis s [' *» 1 d z 



I 

J cos C V ) s I ri""~ C y ) d ': f bacaii ";' 3 ! 

1 1 ) t r a n s f o rm a 1 1 f n i ri! c; 1 1 c* n s 

sf kE de I: a rm I ne I f t h s r* e s us 1' 1 1 n i: e,;=; ^^ 

*-% \ f"}' if' /■,,« ff^-f, , , if ifi "'"'I'' ,■':■ »L i|i|, 

/ Si i^'(f)i(l<'sec^ C;.')' ^3t::e3 



! . 



f,, ,,| ^,, ^, cJJT'i '' ,f ^^ ,,, ,„ ,, Ji _, 

ff 'f Hi'lf^ *I||, ft li 

/' 

I "U (» # -■=:, ,,f:, ;•, .,.>,|i , i ,, ■* ■;:, '!"i ,'l . , i! . I :i. 

" ,j "" 

of c^ne of tfis terjiis sqijsl s a ecjirss >: c:r?it: 

tsrrs'JS I 'li 3:1:;' ise prci-d* (ic:"l:,-j s>n p 

a t a Iji 1 € ,: • T li I s i ss a s c: 1 ;> s c :; t ■; s t 'r=i; f> 

t a b I e ot ■ 1 3 1: eg s- a ! s ., i *:; s fitiiu 1 d " b c: 

iri h li t so e \? e r I s ?' e .3 1 1 y ^ 1 • * c c * s s a r «■ „- " A 

I: ii e t: €3 r r 1 n n cle t' c ci'ri s I d ( u r a ii:; ii q i ^ f ; I 

that csne tlis des*l' vat I ¥e I s iygl'cr!!! > 

wh I '5 hev'fjr I s sic^ii iisi^ori s Z as's t .,, 












/ "" "" '" ' "" 












Wi"!! Set's csf' 









f\ If" C-JE"; '':'l*' If ffl'lS^J H'lrv'. ffi'* .1:»». /■:ili!\ f! ■» ■£? .'PV, *«% , p, ft' m* w.«l /' « '" ™ ?! «... j" , ft fPI»'» * 

li t k» ii a a li'g S£) I J fes I; ii i,s I ; i: 3 ii ii g 1 5 J § 3 t; 3 S K r, I: S C I' ii n f •'1 ^ • 






algof Ithias \'l ' I si^c. IK, Vv l: I' ■■,:;; ■■■...• ,:::>r"v:,: ; 

r ou I: ! n e c s I ll e c i* !F'; ;''" f I'i i .,.::! I ;;;' t -- ■ ■■; -'j ; • .„ :• ■ :; . 

t lii e i* QT'm ! f ■! I n o; : i - 1 ;:;;c.^'^ i- ""i *: I! ■> .■•;:: :;•:; :\:\ ;:: i . :;: ; ; ■: : ' ': ; ,:' , :„ 

1 9 rm s a n cl I oo I' £ f c i s 'c- ■ "'; ':.; I is \ :,:- f. . ". b o >i , ■ ;; ' 

con t a I n I rjg I: r I go :'v :::r;: » " : ;„ • I c fi 5 r.. c. t " : v. ;i • "': ■ ., 

n c s 1 1: li! 3 s f q i ?: ii cl c • ; =! ''"! „, i:,: h e .-- 1' . * :* ;:■ '■ ■• • ? :' ■ ' ■* ! c, •? : : ••• ■ 

T !"i I s r o u t ! s i e f ::, ■ r c: s •; :■ c i"i s; ! '3 1 ;:■; i o ; ■' ' :■ \ • •:; ■: ;. "■ ■- ; 

bi; FORivi I:! ; ecci v ^'iSts. I f ;::>c ^ t: i-r- ;:':': :":*'■:■;- : 

pJ5*ot:)]c!rvi I s i i ::;^:2 •:;:;! :'''*;:£:■ Ij I i' :;;■.:::;;,,, ;^:"' 

II ii. f o J! I: ! ! n a 1: 1 i? f ^^ 3 ,:: I ^ t: f ■ 1, :;:? .:;.■ =; e i' ■ ■ ■ : ■ : , ^ ',. - '' : , 

1 ri I: e g r a t f n p « * c g s -13 ■;.;■; a ;: :. n t r 1 1.: "I .•■ ■■ " ,::: ■ i : •: r„ 

msi de ,1 M Off t h e 1 e £ :s ,p F R: 'v' <;• :;;; (:; t I ? ::;: ; ; t \ i ': ; ■ ; : 

c- :> d 2 , VI e b e i I e v e t fi "z e 3 „; ' , i: i:"',; r (■•; •:" . ; ,.; : c ' ' , ; • ?; j; 

a p r e 1 * r o c *» 3 !~ I r*; g ;■■ ;: ^ {^ 1: I r\ c ■ , i: c '; ; ;■'.;: ■ : ¥ i : ■ : , ": ■ : „.. 



:.; '":■ i i 



The !''?3,id';r' Gmr'av': q iyrc^s :*;:';',:{■'•,;:" ' 

1 11 1 e 15 r a 1: 1 o ni y? y '• / ( l •• >;, ' ) " ■ : • • :> ::: .; ■■ „ ' ; - ■ * ,; 

rots tins vji'ii^ci'i ii,::,:;;' }y'h:r.: cysy 

slriC:;fl«x, ?' :; ; m^iiy y^y'^ I cr • ];' , |f,ijf y, 

pro h I em I a s e n t t c ■ i;;: ' t ;;;} i • ( ,„i i;; f .■'■;: : ■: vr : > f i 

t he s ob s 1: 1 1: t 1 1, i o n s: ^- 1 i ■■:;: •: ; / } I ::, ; :,.. ■': >-.:- : 

p r o b I em I s s s n t t c 14 ;••! i i : ■ v ; : - -■■" I 

integral 1/5 „ -f -^y :*>, - „■;'■;,■;; ,. : sC?: '^^ 

so I union ■"■ 1.^3 Zviflii tr€i< J^;':i { ,?: :':; .^ '•■':■ 

- the s a me s o I s.i i: ! ! ::i rii t h y t if .::3 s > b t r. i v 

Sjy rIT [5; H ■•?• "^rm:}<i V:f:i : ::;';i:,v!' ■■ h ■::!■:' 

SCI lot! o^ri ,, bu ':,; a i so C'l'-i is -■ :; ; ^d r;;. :,;• £ :} ' ■■ ■ :; : '• ■■ ■ 

rc5warcfliig,j ,;: .um e s ■:; yt , ''(-yi , :; i •;;;,: ;^ 

solution,,, Qur QTograyi c'l .i I :,: 'i' , 

- B h a r f n g ., i i : t i: s g d a r:* ::; e c j, :; ! :' ■ i y t y ' " , :•', "„ 

of wh I c ti ? p :; G b iS b 1 v' ;,■■-!■ : ■ :.• { i : :::• n y - „ ;;: " , \ . ;:• i' 

an d ccffRTiiisii I c 1 1; 1 i i g v.;? :*' 1 1' * t t h 



.>:j' 



Coiiisldei!'i ti5|: thij ?4i;cy%s 

3 1- is irali Emd fwedlat ■:;:!:::•::'■ , r; 

near ftitoi 'e t^ilfscii ) '^l 1 'i '':Ty::.r.'t 

s t: a nda r ci i: a b :' e of f i;-; f:,; e-i j ;•■■:: ! ;■ ;■ ' ' 

PB P« 5 ccwTSpi! 1 1 c;; *;' 0, t P I: " ■;) ;': (' ' c, ■• •: f ■:■/?, r 

wc-rd File f!!!!! if! ?■ >:he t i' r;:,;:!:'! .:]:.;;::■■;''■;:: 

vastly redf sc€m:'I £;.;■■::■ 

siilirois itfni !S is. .ii/i bv'' :.:::?'■';- '■^'■^ 

Will 11© SAifl r 15, HkAm t-y'lPr C: -•■;if 

a iarge rrjemori^,, '■■;:::" ':'..■::,:' ■ ^. !,■•'':: 
pe r f o rrae d «■■ s \i5^£. "ij ij -,. ,, .: -,. ,: .■ , ;., :: ,, ,, ,,,:. ., ,, . 

speed a fid I;'; ' the rsf'-yn , '' c 



Oi": 



';: V t" 



The trouble V'*f!th the algorithms for Integration Js that 
they usually deal v^lth classes of expressions v/hich sre 
unique In that each meraber of these classes Is integrabl-a. 
Since many (probably most) problems are not Integrable In 
closed form^ there are likely to be few such classes^ cind 
Indeed few have been found. Hence the value of a purely 
algorithmic program hinges on Its ability to handle very 
large^ but not very complicated^ problems* This Is useful 
but, we feel^ not very Interesting nor dffflcuHo 

The truly difficult problems In Integration are the 
ones whose IntegrabI 1 f ty Is really in doubt. In general 
there can be no decision procedure for Integration* This was 
recently proved by Richardson?^ We ha\'e found a simple 
decision procedure for the case R(x)expCP(x) ) where R Is 
rational and P Is a polynomial. It Is rei:}uted that there 
exists a decision procedure for Integrands which are 
algebraic functions (these Include rational functions and 
roots of rational functions). This Is yer^ surprising^ If 
true. One must resort to heuristics In order to be able to 
solve efficiently the problems that remain^ Vie might try to 
differentiate all possible well formed expressions and match 
these with the Integrand* This Is clearly Inefficient^ and 
may take forever^ One heuristic which we have found (v;hlch 
was motivated by discussions with Professor M, Minsky) Is 
the Edge heuristic oientloned above ^ It Is surprising hovj 
well it performs by Itself. With the addition of a few 
algorithms such as factoring It should be able to solve all 
problems which are Integrable by the algorithms and many 
others* 

The rule which governs the use of the heuristic is as 
follows: determine the most complicated subexpression In 
the Integrand^ guess at a form of the 'ntejgri^l vihose 
derivative contains this subexpression^ then determine the 
value of the variables In the form by making further 
guesses, ^Therc Is theoretical evidence In the ^^lorks " of 
Llouville that this is a very reasonable approach, Ws shall 
describe the guessing rules (as far as we know them -rt 
present) and give some examples below^ 

Below a^b shall stand for undetermined functions of x. 
Tne expression Involving gCx) Is the most complicated 
s.jbeKpresslon In the Integrand^ and fCx) Is the rest' of the 
Integrando We assume that the Integrand Is a product of 
t^srmso 

Form of the Integrand form of the Integral 
f(;«)(gCx))^ 



f CK)exp(g(K)> 



a(g(x)) > b 




n^O 


a(g(x)) ^ 


b 


n< 


aCexp(gCx))) 


4. b 
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fCx)sin(gCx)) acosC g(K)) -J- b 

f(x)cos(gCK)) asinCgCr.)) -s- b 

f(K)logCgCx)) alogCgCx)) * b 

When these forms fall to yield an answar^ then the 
Integral may contain logarithmic terms (I.e.. constants times 
functions of log^ arctan^ or arcsfn)„ 

f(x)gCx) alogCb'i^cgCx)} -i- d 

One should switch to the arctan or arcsin formulations 
when encountering complex constants In the above. 

An example of the use of th^ Edge heuristic* 



X 



r4 



Try a(l « ;?P ^ b 
Differentiation yields 

-3/2 ail - js-^)^(-2k) ^^' a'' C 1 - r.Y^ « b' « x^'ci - / )^ 

Try a « xVC3x) =* 1/3 x^^ hance a '' »' k"^ 

and b' -"« -x'^Cl - x"" )'^ 

in order to solve for b try 
cCl • x')'^ * d 
-1/2 cCl - x-')^C-2k) •«• c'ci « x)^ '^d' « -/C'l -• i}^ 

Try c «'' '*K^/n « -Jt 

/ -4- 

Hence d « CI - x"') 

The solistlon for d Is arcslnCx}^ which Is obtained from 
thm logarithmic term f!ogC-fx •«• CI « x:'^ ri. 



The final ans^^er Is 

1/3 x'ci - K^r-xCl - K^}^ '>arcs\niK} 



We must explain some of the steps above which are not 
obvlousc.^; In this problem the ^hardest sube^tpresslon*' Is 
CI • x-^) -'a This 'hardest subexpression® Is chosen from 
terms in the Integrand (assuming that the Integrand Is a 
product) which can not appear In the derlvatlv^a of the other 
terms* If a logarithm appears onlv oncG In the Integrand^ 
th.Q term containing ft Is a good candidate for the most 



compUcateci subestpresslon* In Zha df vferent latGcl e:cprss3lon 
a or a'' Is selected such thet the terra hwolving ft vdill 
equal the Integrand, We v^ll i choose b^ so that It equals the 
rest of the e;?presslon« 

The two problems which SAINT was able to solva arid 
which the algorithmic program vias not able to solve arc 
xcosCx) and xexp(x)/CI^xr * Both of these problems can bs 
solved with the aid of the Edge heuristic. The solution of 
the latter problem led us to e^ decision pror^ore for th-j 
case R(K)eKp(P(x)) . The prable?Ts cos(*> k) vihSch was 
attempted and not solved by SAIMT can likewise be solvt^d b7 
the Edge heuristic* 

There are two criteria by v/hlch one detarmines when^ to 
terminate the use of a guess made by the Edge heuristic. 
Vlhen one has generated the same subgoal as the original 
problem^ then clearly one must quit the current guess, \^hen 

a subproblem Is generated which Is a constant multiple 
(different from one) of the original problertj, theii 

transposition Is used to obtsin the ariswsrc Vlhen the 
subproblems generated by the heuristic tend to grovv, ons 

should give up the current guess^ altho-jgh this Is a 
situation which Is not so clear ciii\ 

There are son-?e cases with whsch the h'aurlstic w51l 
encounter difficulties* in the case of rational futictfons^ 
the denominator must be factored before the heisrlstlc has a 
chance of succeeding. Otherwise?/ there Is no simple v.;ay to 
make a guess at the 'most complicated sybe>ipresslon® , A 
similar difficulty occurs in thvj case of rational functions 
of trlgonoiYsetrlc functions. The algorithms clearly can ts-ke 
over In these casesc 

Another source of dlfflcyTiy occurs In those problems 
vtfhose solution contains more than one; logarrthmlc terrne Vie 
knovj of no simple way to break up the problem Into 
subproblems which yield the logarithmic terms a This Is tfie 
reason why the. case of algebraic functions Is considered so 
difficult.. 

The rear source of difficulty to an integratlcm 
program^* and to the Edge heuristic In particular^. Is dije to 
the fact that the character of the Integral may not be v/hst 
It seems because some transformation has tak'an place vdifch 
masks ?to Suppose sCx) Is equivalent to zero^ and we are 
asked to integrate itc Unless wa realize the eqislvalence<^ v;e 
may have great difficulties In obtaining; the resolto This 
appears to be thiB real reason vihy Integration ls# In 
general^ recursively unsolvable. Every express! on which !s 
Integrable Is equivalent to an expression which Is trivially 
lntegrable# namely the unslmpllfled derivative of Its 
Integrals The difficulty with Integration can be said to be 
due to the transformation which has taken place In the 
derfvatlvee A single program ca';i only counteract a finite 
number of these transformations., but their number^ In some 
sense^ Is Infinite,-. 
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Finally^ the appUcability of the Edj^e he.jriGtl: ::-.tops 
with linear^ first order^ differential eqaatlc^s. Sf'ice" thi* 
form of the solution of nonllneair differentia] e^nustioms 
critically depends on constants in the eqjatsoji, othf-'f 
heuristics will be necessary to obtain tli« sol-jtlon to rhe'^p 
problems* 



PASE 11 

We shall now describe tvjo of tha basic roti tines ii<^acI by 
the aigorfthmlc Integration program, 

The matching program - SCHATCHEN 

There Is a set of routines In SAUfT v^hSch SIsgle did 
not describe at all. These rout!n^s comprise a rerw-^^ve 
matching program ccslled El Inst CELementarr INSTance) "'" In 
about two pages of LiSP code Slagle wrote a very powerful 
I«Ii !2L f ^^ ^2 recursive that In months of examining 
and using It we have not been able to figure It out 
completely* It Is so powerful that In the six years since' It 
was written no previously documenttKi attempt at a mafchtp'* 
routine has corns close to equaling Its power, r not Formul'a 
^!f?!! Tc^^^'^^r f.1^ simplIfScatloii progra^^ nor the FAMOUS 
system of FenlchelV Most unfortunately, Slagle badly misused 
^vi!. "^ *''®* "® v^as faced with enoniwus difficulties In 
rittJng his program In core and thus many of his subroi^^frrs 
made use of the power of El Inst In order to F^aln space at 
the expense of time. Probably a gaVa of^'an ord^^r of 
magnitude In speed could have been attained by avoldInF 
EHnst In some matching sltuatfons and by usfn^^ sofx-lal 
purpose matches as we have done, .^.>-tja8 

-^vs s-*^\^!^® ^^^^5®^ ^ matching program called SCHATCHFf^ 
'lYlddfsh tor match-maker) which Is a take«of'f* on Ei!nn«- it 

LuZ"}eclrlZZ'''' "' "'''''' '"''■' '' ^-"-^^^"'-^ '"^^ -^ 

SOHATCHEf^ Is a f?inctlon o? tv;o argw^icmts - en 
exprpslon and a pattern. its surpose ! «- ^o ' d^f-rirSn- 
whether or not there exist values for the variables. ^ In '^^^ 
pattern tor which tfm petteai becomes equivalent to rh^ 
ejipresslon, its value Is either MIL or a dfctlosiary n^ f'h- 
'/hl"?%f ^'^^ variables. It Is very similar In purpose to 
the left-hana side of 3 COMiT rule, but Is spr^cl f Icall? 
designed for algebraic eitpressf ens, p-c.i.6cji«/ 

SCHATCHE^J assumes th&t the operators PLUS and TIHFS pr^ 
comntutatlve operators with variable number of argonK^nts'^ 'ft 

ImfX^'Z/Llr' "trf ^f^^^^P^^ Involving and 1 for PLUS, 
HMp and EXPl,: ft Is this fad: which g!ves It and El Inst 

Ui'Othered by massing operators. 

Below we present a somewhat f Ictlo^nallsed (Vsrriotlo'i 
of the program and the patterns that It accepts. -^'-^^-^^^-^^^'^ 

succeldsc ^^'' ^''^''^^''^'''' ®^^^^^^^ ^^'^^ pattern, the match 

If th© pafctern Is of the fom (VAf! a g sr^l ^^rr? '« 
c^^id the e^presslo^ Is e. thsr^ gCe arsl arg2 ^i: V" ,; 
evaluatedo if the ^al^e of g Is 'false, the matih fallC 
Otherw se . the match succeeds aod CCa . e)) fs aopend^d t; 
a dictionary or valsjes,. .^-ui.„u *.o 
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If the pattern is of the form (op pi p2 ««. Pjm) where 
op Is either PLUS or TIMES, then the following takes place* 
If the expression e Is not of the form (op el e2 •«« ak> 
then It Is converted to (op e) ^ 

For each subpattern p^ a search Is mad© of the 
expression for a subexpression ej which will match It* If 
one Is founds the match continues v/Ith the next subpattern^ 
If none Is founds, a match Is attempted between p^ and If 
op*PLliS or 1 If op»TIMES. If the attempt succeeds the match 
continues with the next subpattern* Otherwise the match 
falls. 

If after all the p^ have been matched, the expression 
still contains some terms, an attempt will be made to match 
q with them. Hence q serves the function of the $ In a CO?^IT 
rule. If no terms remain, then q must match If op*»PlUS or 
1 If op-TIMES. 

If op«PLUS and If one of the 9^ Is of the form (TINES* 
rl r2 .«« r^ (VAR a s args)) then what Is desired Is that a 
should become the coefficient of (TIMES ri r2 *,. r;^ ) In 
the sum. This Is done by looping over the remainder of the 
expression and matching (TIMES rl r2 „.. r^ (VAR a s args)) 
with each summand In It. The value of the match Is the sum 
of the Individual matches. 

If the pattern Is of the form (EXPT pi p2) and the 
expression e Is of the form (EXPT el e2), ther, pi must match 
el and p2 must match e2, or pi must match e and p2 must 
match le If e Is 1, then p2 must match or pi must match 1, 
If e Is 0, then pi must match Oe 

If the pattern Is of the form (op pi p2 .«. pk)^ and 
the expression Is of the form Cop el e2 *, *. r ek)^ then each p^ 
must match e j « 

All other matches fallc 

Suppose we wish to match for a linear ejipresslon In Kc 
The following pattern may be usedc 

CPLUSCTIMES* X (VAR B FREE)) (VAR A FREE)) 

Her® FREE Is the name of a routine which checks to see If 
Its argument contains an x. This pattern when matched v^lth 
the expressions belm^ at the left will yield the results at 
th@ rights 

5 CCA . 3){E „ 0)) 

X CCA . 0)CB . 1)) 

(PLUS X PI 2 (TIMES y x)) CCA (PLUS PI 2))(B (PLUS 1 y))) 
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There Is another facility In SCHATCHEN which allows the 

pattern to specify a loop over the expression. Suppose^ we 

were Interested In performing trigonometric simplification 
such as 

a sin (b) '^ a cos (b) •«• e « a * e 

The following pattern will perform the matching 
necessary for the application of this rule, TRUE l<j s 
function which will accept any Input. 

^?iV5St?2?lI!^^^{^^'*^^^5'^^V^« S TRUE))2)(VAR A TRUE)) 
(TIMES(EXPT(COS(VAR C SCHATCHEH 8))2) 
(VAR D SCHATCHEN A)))(VAR E TRUE)) 

With the pattern above a loop will take place over all 
sm and cos-* expressions until one pair Is found which 
meets the criteria or until all are exhausted. We do not 
necessarily recommend this approach to the sfmpUf IcatSop of 
trigonometric expressions* 

It would be useful to make some of the folloivlns- 
e;ttenslons to SCHATCHEN. Currently SCHATCHEN can report a 
succesful match without yielding a dictionary contalnln'^ a 
value for each variable In tfie pattern. This could be 
changed^ for eKamf>le^ by giving defiiult values to the 
variables. On a more semantic level, om? should consider th^ 
VAR pattern to be only of a number of possible mode% of 
f!?5i"^^LSL^''l?^^®^v"**^® ^^^ ^^« corresponds to the BUV 

T^^L''^,,?2'*^^r- P2^^ ^^^"^^ ^^^ "^^^ corresponding 
to the UAR mode of CONVERT. This would accept any value 

which meets some criterion the first time the variable Is 
mSJfr^r^^i ^if* ^^^^ successive value of the variable must 
match the first one. This Is similar to the nym*-r?cal 
const tuent of COMIT which must match a previous 
constituents From an algebraic standpoint SCHATCHEM Is 
limited In that It performs no divlslonc For fnstanr©^ |r 
cannot match a pattern which would correspond to ax *'to an 

?n?Ii«2n!2!] ^"H"^ ff ^ .Some sort of division should he 
Introduced, but Its applicability should be governed by th^ 
user through a new operator such as TIMES**. M wotjfd ^l^-o 
be wise to maintain the mode declarations such as "'vAR 
^rSlicMFM !^® t^'^^^fr"' ^^^^^^ ^^ separate arguments to 
f ^f!f !? ^^ '". ^1^"^^ ""^ CONVERT, or as some global 
declarations as n Formula Algol. Furthermore there should 
e.Klst some facility for performing the construction of an 
expression through the use of the dictionary supplied by the 
matcho This can range from a simple construction routlnr^ m 
a very powerful one as In CONVERTo ^^''^^^ *.- 
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The simplification program - SCHVUOS 

We have written yet another sffapli f Icatlon program^ it 
Is called SCHVUOS which Is an acronysTi for SCHatchen*s 
Version of an Unassuming Operational SImpllfler, It Is 
remarkable In that It Is very short. The basic 
simplification routines are only two pages of LISP code 
long* It assumes no standard form of the expression,. Hence 
It Is Inefficient on ^ery large expressions* However the 
size of an expression occurlng In Integration Is relatively 
smatl# and hence It Is quite fast for our purposes. 

In a sum the summands are simplified^ A summand Is 
first stripped of Its constants and a match Is made for a 
constant multiple of this summand using the Tir^ES* facility 
of SCHATCHEW. This Is done until each summand has been 
accounted for. In a productj» after the sSmpllf Icatlor? of 
all the terms^ the program collects exponents by using 
SCHATCHEN and matching for a variable power of the terms In 
the product. The usual transformations Involving and 1 
with the operators PLOS^ TIMES, and EXPT are all made. 
Furthermore an exponential term whose base Is a product Is 
expanded^ and constants are multiplied through sums. 

Functions such as log and sin have their own 
simplification routines which perform noti-cojutroverslal 
simplifications such as slnCO)«0. 

There Is an Indicator which tells SCHVUOS whether or 
not to simplify subexpressions*, ft Is usually on except 
during differentiation where the differentiated expression 
Is built from the bottom up and Is simplified at each level, 
with no redundant simplifications being performed. /o^Th!s Is 
SCHVUOS* alternative to the AUTSIM bit of FORMAc"'^ and to 
Martin's version thereof^v 
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The following routine Is not part of the current 
algorithmic Integration program* 

An Integral TAble Look-Up - I TAIL! (read l-te11-*yoy) 

We had spent some time building an Integral table 
look-up before realizing that for Integration It would quite 
likely be useless since the program could^ In all 
probability^ Integrate anything In the table quickly. Since 
this Is not the case for differential equatlonsj, definite 
Integration^ or summation of series^ we Vtould like to 
present what we have done In ITALU* 

The basic steps In ITALU are as follows: An express ton 
to be Integrated Is first hash-coded. The code (a floating 
point number) Is looked uq In an array (binary search) and a 
disk address for It Is found. The disk Is read and forms of 
the Integrands In It are matched with the given expression 
(probably only one exists with the right code) until one Is 
found which matches the expression* The coresponding 
Integral Is evaluated after substituting In It the results 
of the match. 

For example^ the hash-code for CEXPT £ (PLUS Pi (TIMES 
2 x))X)) v/ould be the same as for the form 

CEXPT(VAR C FREE)(PLUS(TIMIiS*-(VASl ¥ VARPXVAR 8 FREE)) 
(VAR A FREE)) 

(where VARP tests to see whther Its argument Is equal to x) 
and the Integral found by ITALU would then be 1/2 CGXPT E 

(PLUS 1 PI (TIMES 2 X))). The trick of evaluating the 
integral allows the integral to be a function and thus 
allows Iterated integrals to be easily obtalnedc 

The hash-code Is designed to check for the form of the 
expression and to Ignore constants with respect to the 
variable of Integration. Thus k Is coded like a^br. , but 
s!n(2x) Is coded differently from sfn(x'^) or cosC2k)„ Coding 
Is done recursively* The code of a sum Is the sum of the 
codes of the summands. Ignoring cons tan tSo The coda of a 
product Is likewise the product of the code of the terms^ 
Ignoring constants* Trigonometric snd logarithmic functions 
are coded by exponentiating the code of the argymants by a 
constant which is different for each function* The cods for 
an exponential term Is more Involved^ Since the exponents 
"•1.-1/2^ -1/2^ 2^ -2 occur so frequently^ these are made special 
cases© All other constant exponents are given the same 
codec The code for CEXPT a b) Is the code of a raised to the 
code of be 

We have also experimented with a hash-code which notes 
It 8 constant occurred^ but ignores the value of this 
constant* Thus l^-x codes differently from x^ but the same 
ar» 2'e'Xo 

The advantage of the schema above Is, that It Is quite 
fast (we believe that one can get running times of about ^ 
second per Integral^ most of which Is spent accessing the 
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disk)* Furthermore the speed of the look-up Is not 
appreciably deteriorated with an Increase In the number of 
entries In the table* The scheme yields a powerful table 
look-up because of the use of SCHATCHEN to perform the 
match^ which allows variables to appear In th& Integrand 
just as they do In the standard Integration tables* The use 
of the device of evaluating the Integral allovis Iterated 
Integrals to be conveniently entered In the table* 

Integral table look-up schemes reported In the 
literature are to be found In R 
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